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Abstract. In this paper we present a statistical analysis of the 
7-ray spectra of flat-spectrum radio quasars (FSRQ) compared 
to those of BL Lacs. The average spectra and possible system- 
atic deviations from power-law behaviour are investigated by 
summing up the intensity and the power-law fit statistic for 
both classes of objects. We also compare the time-averaged 
spectrum to that at the time of 7-ray outbursts. 

The spectrum of the average AGN is softer than that of 
the extragalactic 7-ray background. It may be that BL Lacs, 
which on average have a harder spectrum than FSRQs, make 
up the bulk of the extragalactic background. 

We also find apparent cut-offs at both low and high ener- 
gies in the spectra of FSRQs at the time of 7-ray outbursts. 
While the cut-off at high energies may have something to do 
with opacity, the cut-off at low energies may be taken as indica- 
tion that the 7-ray emission of FSRQs is not a one component 
spectrum. 

Key words: Methods: statistical; BL Lacertae objects: gen- 
eral; quasars: general; Gamma rays: observations 



1. Introduction 

Up to the end of 1994 more than 50 extragalactic radio sources 
were detected as emitters of high-energy 7-rays by EGRET. 
The majority of the sources are flat-spectrum radio quasars 
(FSRQ) and about 20% are classified as BL Lacertae (BL Lac) 
objects. Especially in their flare state FSRQs often emit the 
bulk of their luminosity in the form of 7-rays (von Montigny 
et al. 1995). For a number of objects COMPTEL data at a few 
MeV (Schonfelder et al., 1996), OSSE data at a few hundred 
keV (McNaron- Brown et al. 1995) or the extrapolation of X-ray 
data imply a spectral break at a few MeV photon energy. This 
break may be attributed to incomplete cooling of radiating 
electrons (Dermer and Schlickeiser 1993, Sikora et al. 1994) 
or bremsstrahlung and annihilation emission of cooling pair 
plasma (Bottcher and Schlickeiser 1995). 

Individual 7-ray spectra in the EGRET range, i.e. above 
30 MeV, can generally be well described by power-laws. For 
the BL Lac Mrk421 the power-law behaviour extends up to 
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TeV energies (Punch et al. 1992). Though a number of pa- 
pers deal with the multifrequency spectrum of 7-ray AGN, 
to our knowledge no attempt has yet been made to investi- 
gate whether there are systematic deviations from power-law 
behaviour in the 7-ray spectra of AGN. Due to the different 
environmental conditions one might expect a different impact 
of opacity effects, for example, on FSRQs on one side and on 
BL Lacs on the other side. 

In this paper we analyse the class-averaged spectra of AGN 
by summing the observed intensity and the statistic of power- 
law fits to the observed emission. We also derive the spec- 
trum of the average 7-ray AGN which is to be compared to 
the spectrum of the diffuse extragalactic background. In the 
next section (Sec. 2) we consider the time-averaged behaviour 
of sources by using the summed EGRET data of Phases 1,2 
and 3 corresponding to observation times between May 1991 
and October 1994. Since many of the extragalactic sources are 
variable at 7-rays we discuss the behaviour at flare state in 
Sec. 3. A discussion of the results is in Sec. 4. 



2. The average source 

Here we have summed all EGRET data taken between May 
1991 and October 1994 to obtain 7-ray allsky maps of counts 
and exposure. A quality cut has been imposed by restricting 
the field-of-view in the individual viewing periods to 30° off- 
axis angles before summing. At energies above 100 MeV we 
have performed a standard search for point sources using a 
maximum-likelihood technique (Mattox et al., 1996). The data 
in the energy range between 30 MeV and 100 MeV are usu- 
ally neglected in this task since they add little to the statis- 
tics and source positioning, but require the energy-dependent 
point-spread function be averaged over an even larger en- 
ergy range with corresponding systematic uncertainties. In 
total we found 148 excesses with a likelihood test statistic 
TS — (2 In Ao — 2 In A) of least 9, of which 44 are positionally 
coincident with AGN known to emit 7-rays in the EGRET 
range. Table 1 lists these sources which can be further divided 
into BL Lacs on one side (11 objects) and FSRQs on the other 
side (33 objects). All sources except one (2155-304) have been 
observed with likelihood test statistic higher than 25 in individ- 
ual viewing periods or specific combinations thereof (Thomp- 
son et al. 1995, 1996). The source 2155-304 is observed with 
likelihood test statistic of 33 in a later viewing period (VP 
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404) during Phase 4 of the CGRO observing program (Ves- 
trand et al. 1996). We do not intend to find new AGN in our 
sample, but select the already identified ones for further anal- 
ysis. Our threshold criterion of TS > 9 thus corresponds to a 
detection significance of at least 3 a in the time-averaged data. 
Some AGN which have been found with high significance in in- 
dividual viewing periods fail to show up in Table 1 since they 
have been absent in other viewing periods so that their average 
signal is below threshold. 
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Fig. 1. The summed intensity spectrum of the 33 identified FSRQ 
and 11 identified BL Lacs in the total EGRET allsky data above 
100 MeV. The error bars are derived by Gaussian error propagation 
of the individual uncertainty measures. The spectrum of FSRQ is 
significantly softer than that of the diffuse background, while that 
of BL Lacs is consistent with it. 

For all 44 AGN we have performed a spectral analysis. The 
positions of the 7-ray sources have been set to those of the cor- 
responding radio sources before doing a likelihood analysis of 
all 148 excesses in ten energy bands. Since the best model for 
the data consists of diffuse emission and 148 sources, we have 
to fit all these quasisimultaneously although we are interested 
only in the results for the 44 AGN. Observations of 7-ray pul- 
sars have shown that the calibrated effective area in the two 
low energy bins (30-50 MeV and 50-70 MeV) is overestimated. 
In the standard analysis the observed source flux is multiplied 
by a correction factor which accounts for the miscalibration. 
The accuracy of this correction factor can be estimated to be 
better than 30% at 30 MeV to 50 MeV and better than 10% 
at 50 MeV to 70 MeV (Fierro 1995, §3.4). 

We have summed the observed intensity in the ten energy 
bands to derive the spectrum of the average AGN. This spec- 
trum is what we would get as contribution to the diffuse ex- 
tragalactic background if the AGN were unresolved. This 'av- 
erage' AGN spectrum can not be directly compared to other 
studies in which power-law fits of the spectra of individual 
AGN are averaged (e.g. Mukherjee et al. 1997). However for 
the purpose of comparison with the diffuse extragalactic back- 
ground co-adding of the observed intensity spectra without any 
normalization or fitting procedure is the appropriate method. 



Table 1. A list of radio sources which can be identified with 
point sources in the summed EGRET data of the time pe- 
riod 5/91 to 10/94 at energies above 100 MeV. The sources 
are categorized to belong either to the BL Lac class or to the 
class of FSRQ. The second column indicates the EGRET ref- 
erence with acronyms for the second EGRET catalogue (2EG, 
Thompson et al. 1995), its supplement (2EGS, Thompson et 
al. 1996), and the papers by Mukherjee et al. (1997) (M) and 
Vestrand et al. (1996) (V). If there is at least a moderate level 
of variability, the third column gives the viewing period num- 
ber for the peak flux in standard EGRET notation. No entry 
in this column indicates no evidence for variability. Note the 
following constraints for the selection of the viewing period 
of peak flux: the total statistical significance is required to be 
> 4 a and preference is given to data taken with smaller aspect 
angle in case of viewing periods of comparable duration and 
comparable source flux. 
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Since the radio catalogues are incomplete near the galactic 
plane we do not expect - and in fact we do not get - identifica- 
tions for sources located at |6| < 10°. Thus effectively there are 
only ~10.4 steradian sky area which we have searched for point 
sources and which is to be used to derive the corresponding av- 
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erage 7-ray intensity. The results is shown in Fig.l separately 
for FSRQs and BL Lacs. 

The sum of the 7-ray intensity contribution of all 44 AGN 
is about 7% of the total diffuse background intensity. It is in- 
teresting to see that the 7-ray spectrum of the average BL Lac 
seems to be harder than that of the average FSRQ. The differ- 
ence in spectral index is Ss — 0.12 ± 0.08 and is thus not very 
significant which is mainly due to the large uncertainties in 
the average BL Lac spectrum. This does not imply that in sin- 
gle viewing periods BL Lacs have always harder spectra than 
FSRQ. In fact we see a remarkable spread of spectral indices 
for both classes of objects when individual viewing periods are 
considered (Mukherjee et al., 1997). We also know that indi- 
vidual sources can change their spectrum: there is a trend that 
the spectrum hardens with increasing flux level (Miicke et al. 
1996). But this concerns individual sources. The spectrum of 
the average in both object classes appears to be different, and 
therefore they would contribute with different spectral charac- 
teristic to the diffuse 7-ray background. 

The spectrum of the average AGN is dominated by that of 
the FSRQs and it differs significantly from that of the observed 
diffuse cxtragalactic background (Kniffen et al. 1996; Sreeku- 
mar et al., 1997) as we show in Fig. 2. When fitting the ratio of 
both spectra with a power-law we find the background spec- 
trum harder than that of the average AGN by 5s — 0.15 ±0.04 
when all energy bins are considered and by Ss = 0.17 ± 0.045 
when only the data with very good statistics between 70 MeV 
and 4 GeV are considered. The goodness-of-fit for a constant is 
8- 10 -3 and 3.5- 10~ 3 , respectively, which corresponds to about 
2.7a significance. A Fischer-Snedecor-test indicates that with 
3.5<r significance a linear relation is a better fit to the intensity 
ratio than a constant. Especially between 70 MeV and 4 GeV 
the uncertainty of the background intensity is so small that the 
assumption of Gaussian statistics for the uncertainty of the in- 
tensity ratio can be regarded as a reasonable approximation. 
Since we see some objects only at flare state when they tend to 
have harder spectra, the true average spectrum of AGN may 
be even softer and thus further away from agreement with the 
observed spectrum of the 7-ray background. 

There is no cut-off visible in the 7-ray spectra with the 
possible exception of a weak deficiency below 100 MeV for the 
FSRQs which may be the outer extension of the usual roll-over 
at a few MeV (Schonfelder et al. 1996). However, the expected 
spectral form of the average intensity is not necessarily power- 
law. Here we have summed over many sources with different 
spectra. Also some of the individual sources have changed their 
7-ray spectrum over the last years. A summation over differ- 
ent power-law like spectra results in a positive curvature of 
the average spectrum (Brecher and Burbidge 1972). Given the 
spread in spectral indices for the 44 sources we can estimate 
that the averaging should change the spectral index by around 
0.1 between 100 MeV and 10 GeV. 

We have also searched for systematic deviations from 
power-law behaviour in the 7-ray spectra of the 44 AGN. For 
each individual AGN we have fitted a power-law spectrum to 
the data. The difference between this fit and the measured in- 
tensity in the ten energy bands weighted by the observational 
uncertainty, i.e. x = (Ifit — I)/(8I), has been summed for all N 
FSRQ and BL Lacs, respectively, to obtain the net deviation 
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where the applicability of Gaussian statistics has been assumed 
for the renormalization -^=. 
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Fig. 2. The ratio of the average intensity of all observed AGN to 
that of the extragalactic diffuse 7-ray background. This ratio is not 
compatible with a constant and thus the observed AGN would, if 
they were unresolved, give a softer diffuse emission than we observe 
in the background. This implies that the background can not be the 
superposition of unresolved AGN with the same characteristics as 
the observed objects. 
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We have verified with Monte-Carlo simulations that there 
are no significant deviations from a Gaussian distribution of 
the number Xtot- The results are shown in Fig. 3. No significant 
deviations from power-law behaviour in the AGN spectra are 
observed, either for FSRQ or for BL Lacs. The BL Lac fit 
statistic is noisy. 
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Fig. 3. The summed power-law fit statistic of the 33 identified FSRQ 
and 11 BL Lacs in the total EGRET allsky data according to Eq.l. 
There is no significant deviation from power-law behaviour in the 
average FSRQ and BL Lac spectrum. The two extreme values with 
opposite sign between 100 MeV and 300 MeV in the BL Lac result 
are unlikely to be real. 



Fig. 4. The summed intensity spectrum of the 28 variable FSRQ and 
6 variable BL Lacs at the time of the highest flux level above 100 
MeV. The error bars are derived by gaussian error propagation of 
the individual uncertainty measures. The FSRQ spectrum at peak 
is significantly harder than the spectrum of the average FSRQ in 
Fig.l. 



3. The peak spectra 

A large fraction of AGN is variable at 7-ray energies. Any 
cut-off arising from opacity effects will be more prominent at 
high flux levels since then the intrinsic photon density of the 
source is high. We have therefore chosen a subsample of AGN 
for which at least a moderate level of variability can be found. 
These sources are indicated in Table 1 by the remark 'yes' 
in the variability column. The light curves of the individual 
sources have been taken from the second EGRET catalogue 
(Thompson et al., 1995) and from the standard viewing period 
analysis of the EGRET group at MPE. In total we are left with 
28 FSRQs and only 6 BL Lacs which are variable. The analysis 
is now similar to that described in the previous section except 
that the spectra are not derived on the basis of the summed 
data of Phases I-III but only on data of the viewing periods 
in which the sources showed the highest flux levels. We should 
note that detections with formal significance < 4 a have not 
been considered. In case of two or more viewing periods of 
comparable duration and comparable source flux preference 
has been given to data taken with smaller aspect angle. The 
results for the summed spectra are presented in Fig. 4 while the 
summed power-law fit statistic is shown in Fig. 5. 



There is not very much change compared to the average 
behaviour in case of BL Lacs. One has to keep in mind that we 
are now left with 6 objects and the statistics are not sufficient 
to distinguish general trends from pathological individuals. 



The behaviour of FSRQs during their peak phase is more 
interesting. At first we see that the flare spectra are harder than 
the average, at least between 100 MeV and a few GeV. This is 
a confirmation of a claim by Miicke et al. (1996) who found a 
hardening of the 7-ray spectra with increasing flux level for 8 
highly variable AGN. When fitting a power-law to the intensity 
ratio of average FSRQ at peak and total average FSRQ we 
find the spectrum of the average FSRQ at peak harder than 
that of the time-average by Ss = 0.19 ± 0.05 when all energy 
bins are considered and by Ss = 0.18 ± 0.06 when only the 
data with very good statistics between 70 MeV and 4 GeV 
are considered. The goodness-of-fit for a constant is 2.5 • 10 -3 
and 7.5 • 10~ 3 , respectively, which corresponds to about 2.7er 
significance. A Fischer-Snedecor-test indicates that with 4. Oct 
significance a linear relation is a better fit to the intensity ratio 
than a constant. 
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Fig. 5. The summed power-law fit statistic of the 28 variable FSRQs 
and 6 variable BL Lacs at the time of the highest flux level. There are 
deviations from power-law behaviour in the peak FSRQ spectrum 
at high energies and at lowest energies. The BL Lacs at peak do not 
show a signature of spectral structure. 

We also see that at energies below 70 MeV and at ener- 
gies above 4 GeV the peak spectra show some evidence of a 
cut-off which is more prominent in summed fit statistic than in 
the summed intensity. It should be noted that the error bars 
for the summed intensity are mainly determined by sources 
with poor statistics, i.e. short exposures or large off-axis an- 
gles. These sources usually contribute very little to the summed 
fit statistic since both at low and high energies the uncertainty 
in counts exceeds the number of detected counts. Thus the \i 
that goes into the summation in Eq.l is practically restricted 
to — 1 <C Xi < 1- This is different for sources with good statis- 
tics for which even at highest energies the number of observed 
and expected counts is roughly between 0.5 and 10. Thus the 
summed fit statistic is dominated by sources with good statis- 
tics and its results can not be compared directly to the summed 
intensity spectra. 

At high energies a cut-off would be the expected behaviour 
since the strength of opacity effects depends directly on the 
flux level. Though the low-energy deficit is significant at 50-70 
MeV the average flux per source at flare state is still a factor 
of 2.0±0.6 higher than in the time-averaged case (cf. Fig.l). 
The flux per source at 30-50 MeV is the same as in the time- 
averaged case. Thus it appears that the 7-ray flux below ~ 50 
MeV does not change, or at least not in phase with the 7-ray 
flux above 100 MeV. 

If the 7-ray flaring spectrum of FSRQ is indeed a power- 
law above 70 MeV, which is then cut off at a few GeV due to 
increasing opacity, we may underestimate the significance level 
of this cut-off in our analysis. The reason is that the power- 



law fit tries to fit both the power-law part and the cut-offs, and 
therefore will underestimate the power-law part at medium en- 
ergies which is in principle the null hypothesis in the statistical 
test for cut-offs. This effect can be observed in Fig. 5 as a strong 
preponderance of positive deviations between 100 MeV and 4 
GeV. 

To get a better idea of the significance level of the cut- 
offs we have repeated the power-law fits for the peak phases 
of FSRQs under the constraint that now the fit is based on 
the energy band of 70 MeV to 4 GeV and then extended to 
calculate the true deviations in the outer energy bands. The 
result is shown in Fig. 6. At energies below 70 MeV there is a 
deficiency of intensity compared to power-law behaviour with 
total statistical significance of 3.0ct while at high energies above 
4 GeV we observe an intensity deficit with 2.6(7 significance. As 
mentioned before the statistical significance has to be inferred 
from the summed fit statistic and it cannot be easily compared 
to the summed intensity spectrum. 
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Fig. 6. The summed power-law fit statistic of the 28 variable FSRQ 
at the time of their peak flux level. In contrast to Fig. 5 here the 
power-law fit is based only on the data between 70 MeV and 4 
GeV, i.e. omitting the outer energy bands for which a deviation was 
suspected. The total statistical significance of the spectral breaks is 
3. Oct at low energies and 2.6(7 at high energies. 

This result is stable against the choice of source. The high- 
est flux level of a source does not necessarily imply the highest 
significance. We have tested whether our analysis is influenced 
by data with poor statistics by excluding all sources with for- 
mal significance of less than 6 a in the viewing period of peak 
flux. There was no change in the summed fit statistic beyond 
small statistical fluctuations. On the other hand, some sources 
have strong secondary maxima in their 7-ray light curves or 
have been observed at similar flux levels in adjacent viewing 
periods. To do the opposite test we have extended our data 
base by including all observations which yielded an integrated 
flux of S > 10 -6 cm -2 sec -1 or for which the derived flux was 
within 2 a of the peak value. Again our conclusion remained 
unchanged. The problem with the latter test is that though 
the statistical basis improves, individual sources, which appear 
three or four times, start to influence the result. The fact that 
our result is indeed robust and does not depend on a specific 
selection criteria, gives us confidence that it is real. 
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3.1. Possible Systematics 

Pulsar data show that at low energies the effective area of 
EGRET needs to be decreased from the value determined by 
the pre-flight calibrations. The accuracy of this correction fac- 
tor can be estimated to be better than 30% at 30 MeV to 50 
MeV and better than 10% at 50 MeV to 70 MeV (Fierro 1995, 
§3.4). There is another source of uncertainty at low energies 
which has to do with spill-over. At 7-ray energies below 100 
MeV EGRET's effective area decreases rapidly (Thompson et 
al. 1993). Hence photons with original energy of around 100 
MeV, which are misidentified as of lower energy, get a strong 
weight and mimic a high intensity at low energies. For the typ- 
ical E~ 2 spectrum around 10% of the low energy photons are 
due to this spill-over. Standard EGRET software for spectral 
fits accounts for this spill-over as far as the relation between 
counts and intensity is concerned. However this effect also influ- 
ences the effective point-spread function at low energies which 
is calculated under the assumption of an E~ 2A spectrum. Since 
the possible error in the point-spread function is around 10% 
we expect corresponding uncertainties in the likelihood fits of 
similar order. In total we estimate less than 50% systemati- 
cal uncertainties in the intensity level at lowest 7-ray energies 
and around 20% in the second energy band which is to be 
compared to the total statistical uncertainty in the summed 
intensity spectrum of quasars in Fig. 4 where we have 130% at 
30 MeV to 50 MeV and 32% at 50 MeV to 70 MeV. In both 
cases the statistical uncertainties are much larger than the sys- 
tematic uncertainties so that the former are a fair measure of 
the total uncertainty at low 7-ray energies. 

At high 7-ray energies the spill-over does not play a strong 
role, since EGRET's effective area changes only weakly with 
energy in this range. Nevertheless possible calibration errors 
of effective area and point-spread function may influence the 
result. Also if a significant fraction of the 7-ray sources is 
misidentified, then by using the positions of the radio sources 
we may underestimate the number of observed photons at high 
energies where the point-spread function is most narrow. Here 
it is instructive to compare Fig. 6 to the corresponding result 
for the average emission of FSRQs in Fig. 3. There we had no 
signal for a high-energy cut-off. We have tested this also for the 
case that the power-law fit is derived only between 70 MeV and 
4 GeV (as in Fig. 6) with similar result: there is not even a ler 
indication for a cut-off beyond 4 GeV. Since any systematic ef- 
fect should influence also the average FSRQ spectrum we can 
be sure that the observed high energy cut-off in the peak FSRQ 
spectrum is real. A misidentification of the 7-ray sources with 
FSRQ in general, such that the 7-ray spectra discussed here 
are not those of FSRQ, has been found extremely unlikely in 
statistical studies (Mattox et al. 1997). 

We have tested the reliability of our method by Monte- 
Carlo simulations. These simulations would detect systematic 
problems in the analysis tools, which may arise from the small 
photon numbers both at low and at high 7-ray energies. We did 
not detect significant systematic deviations from a Gaussian 
distribution of the variable \tot of Eq.l. 

4. Discussion 

4-1. The deficit at low energies 

We have seen in the previous section that in case of FSRQ the 
7-ray spectra deviate from power-law below 70 MeV. This is 



most significant with 3. Oct when the emission of variable sources 
at flare state is considered. In the time-averaged spectrum the 
effect is less pronounced which indicates that if there is a deficit 
in the average spectra, it will be weaker since we have better 
statistics there. 

The extrapolation of X-ray data or the OSSE spectra for 
many sources imply a spectral break at a few MeV energies. 
These breaks have to be extended features regardless of their 
origin, even if an annihilation feature is superimposed. This is 
because the standard radiation processes like inverse Compton 
scattering are not monochromatic, when the target photon dis- 
tribution is neither beamed or monochromatic. So even if there 
were a sharp break in the electron spectrum, the correspond- 
ing break in the 7-ray spectrum would be smeared out. It is, 
however, questionable whether this is sufficient to account for 
the observed deficit below 70 MeV, which is a factor 10 higher 
in energy than the typical break energy. EGRET's statistics 
are limited in this energy range, and hence it needs some ef- 
fort to produce such a strong signal at 50-70 MeV. We would 
then also expect a signal at 70-100 MeV where the statistics 
are much better. 

We therefore prefer to interprete our result in the sense that 
the 7-ray emission of FSRQ between a few hundred keV and 
a few GeV is not a one-component spectrum, but rather the 
superposition of different emission processes. 

If the jet material consists of pair plasma we may expect an- 
nihilation and lepton-lepton bremsstrahlung to dominate the 
spectrum up to a few tens of MeV while inverse-Compton 
scattering is more efficient at higher energies (Bottcher and 
Schlickeiser 1995). The comptonisation part of the spectrum 
is related to the electron spectrum at high particle energies 
where the physical reaction timescales are shortest. This may 
explain why at flare state we see a deficit below 70 MeV. 

In case the jet is made of ordinary matter we have to con- 
sider bremsstrahlung and Coulomb interactions besides the in- 
verse Compton scattering. 

Generally the observed behaviour can also be caused by a 
low-energy cut-off in the electron injection spectrum (Bottcher 
and Schlickeiser 1996). 

4-2. The cut-off at high energies 

BL Lacs are defined as having only weak, if any, optical lines. 
That translates to them having only little backscattering mate- 
rial in the vicinity of the central machine. Hence we can expect 
different environmental conditions in BL Lacs and quasars, 
as far as opacity is concerned. The main source of opacity - 
photon-photon pair production - requires a sufficient number 
of target photons of the energy 



where fi is the collision angle cosine in the observers frame, 
and E 7 is the energy of the 7-ray to be absorpted. 

Standard calculations show that the interaction of 7-rays 
with themselves is reduced due to the equivalent of the Klein- 
Nishina effect, which applies when the target photon energy is 
very much larger than the limit given by Eq.2 (Pohl et al. 1995, 
Dormer and Gehrels 1995). Interactions of 7-rays with the X- 
ray continuum may affect the total 7-ray spectrum. What we 
observe, a cut-off only at highest energies, looks more like an 
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additional target photon source coming into play. This could 
be an accretion disk, accretion disk emission backscattered by 
thermal clouds, or synchrotron photons produced in the jet 
itself. 

For 7-rays with dimensionless energy ei the optical depth 
for pair production by collisions with ambient photons of en- 
ergy e is 



r 77 = 2n J dz J dfi (1 — m) / de an ph 



where 



16 v 



P 2 ) 



(3- 



/3 4 )lni±| 



2/3(2- 



P 2 ) 



(3) 



(4) 



with 



eei(l - fi) 



For the original accretion disk photon the collision angle is 
unfavorable (/i w 1) so that the UV photons can only interact 
with 7-rays of 100 GeV energy or more (Bottcher and Dermer 
1995). In the direct comptonization of the accretion disk pho- 
tons the Klein-Nishina cut-off may lead to spectral turn-overs 
at around 10 GeV (Bottcher et al. in prep.): that is a basic 
characteristic of the emission process and can not account for 
cut-offs which occur only at flare state. 

In case of the backscattered accretion disk photons 



n v h 



6(4tt) 



■*(e-e ) 



(5) 



a target photon energy e « 50 eV is sufficient to produce a cut- 
off at a few GeV observed 7-ray energy. Here the delta-function 
is a simple approximation of the accretion disk spectrum with 
luminosity L and rjz c / z 2 <JC 1 is the assumed re-scattering rate 
with scattering optical depth r rs = r/z c /z. Then 



Z/46 T r . 
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V e o e i 
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(6) 



with z in units of 0.01 pc and r rs ,-5 in units of 10~ 5 . If T rs ,-5 
is larger than unity, r 77 increases sharply beyond unity and 
produces a cut-off in the 7-ray spectrum at around 6 GeV 
(before redshift). However, correlations between optical depth 
and 7-ray flux above 100 MeV can only occur when the 7- 
ray outburst is due to an increased level of target photons, for 
instance an accretion disk flare. Any variation in the spectrum 
of radiating electrons in the jet will have no impact on the 
optical depth. 

If electrons with Lorentz factors of ~ 10 6 do exist in the jets 
(as is required for the TeV emission of BL Lacs), then these will 
produce a synchrotron continuum up to X-ray energies, which 
provides a large number of target photons for pair production. 
The optical depth for this process is 



L47 



V 5 GeV J 



(7) 



where S is the Doppler factor of the jet in units of ten, rid 
is the effective jet radius in light days (3 • 10 lj cm) in the jet 
frame, and L47 is the vF(y) synchrotron luminosity at energy 
2S 2 (m e c 2 ) 2 /E~{ in units of 10 47 ergs/sec. So if a 7-ray outburst 
is caused by the injection of a large number of high-energy elec- 
trons, the 7-ray emission above a few GeV may be self-damped. 
However, detailed simulations of cooling pair plasma includ- 
ing first order and second order Fokker-Planck coefficients for 
all interaction processes (Bottcher, Pohl, and Schlickeiser, in 
prep.) show that the high end of the synchrotron spectrum is 
likely to be swamped by synchrotron-self-Compton emission, 
which may have a hard spectrum in the relevant energy range 
so that L47 in Eq.7 increases with e, i.e. decreases with E 1 . 
In that case the cut-off would be very smooth and the energy 
at which it occurs would vary strongly for different X-ray lu- 
minosities and Doppler factor. Thus there would be no good 
argument why the bulk of FSRQ considered here should show 
the cut-off at roughly the same energy. 

To summarize: though opacity can in principle account for 
cut-offs at a few GeV photon energy, and seems to be required 
to explain the correlation with the 7-ray flux level, there is 
no clear answer as to the source of the target photons for the 
photon-photon pair production. However, our findings will fur- 
ther constrain models and simulations for the evolution of pair 
jets. 



4-3. The relation to the diffuse •y-my background 

We have seen that the time and source averaged 7-ray spec- 
trum of AGN is softer than that of the extragalactic 7-ray 
background. The spectrum of the average AGN is dominated 
by that of FSRQ, partly since we have three times more ob- 
jects of this class than BL Lacs, partly since individual FSRQ 
are on average brighter. It is interesting to see that the aver- 
age BL Lac has a spectrum which is compatible with that of 
the 7-ray background and harder than that of FSRQ. If un- 
resolved AGN would indeed be responsible for the bulk of the 
7-ray background, the BL Lacs may have to play a stronger 
role than previously thought. Here we outline a scenario which 
would account for this. 

The BL Lacs have on average a much smaller redshift with 
values between 0.031 and 0.94, while more than 50 % of the 
objects in the FSRQ class have redshifts in excess of 1.0. This 
indicates that in case of FSRQ we observe a fair range of the 
luminosity function directly, in contrast to the BL Lac case 
where we see only the tip of the iceberg. In other words, we 
expect the 7-ray log N/ log S distribution of BL Lacs to peak at 
lower 7-ray fluxes than that of FSRQ. As a result the contribu- 
tion of BL Lacs to the diffuse extragalactic 7-ray background 
may be strong despite the small number of directly observed 
objects, and hence it may be that BL Lacs provide the bulk of 
the 7-ray background. This would of course require that FSRQ 
and BL Lacs do not have drastically different evolution prop- 
erties. The reader should note that at least for radio selected 
BL Lacs co-adding of undetected sources has resulted in a ex- 
cess of around 3cr significance (Lin et al. 1997), which provides 
further evidence that this class of sources does emit 7-rays at 
lower flux levels. 

There is yet another effect by which BL Lacs can contribute 
to the extragalatic 7-ray background. Distant BL Lacs may 
emit 7-rays at energies higher than 100 GeV like the close- 
by objects seen by Whipple. These 7-rays will pair produce 
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on the low-energy extragalactic background before reaching 
us. Since the pairs will immediately comptonize the microwave 
background, the energy of the TeV 7-rays will go into an elec- 
tromagnetic cascade and finally reappear in the form of 7-rays 
at only a few GeV energy (Coppi and Aharonian 1997). De- 
pending on the luminosity distribution between GeV and TeV 
emission of the average BL Lac, this process may dominate 
over the direct contribution to the GeV background for partic- 
ular redshifts 2. This effect would results in a hardening of the 
background spectrum in the EGRET energy range. 

We should however not completely exclude the possibil- 
ity that the spectrum of the extragalactic 7-ray background 
is ill-determined. The separation of background and galactic 
emission due to cosmic ray interactions with thermal gas is 
relatively easy to do by correlation between intensity and the 
gas column density. In contrast the expected spectrum of the 
inverse Compton emission is similar to that of the extragalac- 
tic background and the intensity varies little with position at 
high latitudes, so that we we cannot exclude the possibility 
that part of the background intensity at higher 7-ray energies 
is misidentified inverse Compton emission. 
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